Light Curves of Microlensed Type la Supernovae 
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>- ■ A detailed description of the apparent light curves of microlensed SNe la as 

extended and expanding sources is presented. We show that microlensing amplifi- 
\Q • cation can have significant effects on supernova observation. A model light curve 

is used to compare lensed and unlensed cases and we find that significant changes 
in shape can occur because of microlensing. We briefly discuss the probability 
. of observing such effects as well. Throughout the paper we consider source and 

deflector redshifts at 1.0 and 0.05 respectively, and limit our work to spherically 
q ■ symmetric deflectors. 
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1. Introduction 

Out of a number of distance indicators, supernovae have emerged as the most promising 
standard candles. Due to their significant intrinsic brightness and relative ubiquity they can 
be observed in the local and distant universe. Several teams including the High-z Supernova 
Search (Schmidt et al. 1998) and the Supernova Cosmology Project (Perlmutter et al. 1999) 
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have been searching for supernovae at higher redshifts since the early 1990's. Light emitted 
from these 'standard candles' is subject to lensing by intervening objects while traversing the 
large distances involved (Kantowski, Vaughan, & Branch 1995); the further the light source, 
the higher its chance of being significantly lensed. In fact, for cosmologically distant sources, 
the probability is high that a distant point source will be 'imaged' (Press & Gunn 1973; 
Bourassa & Kantowski 1976; Wyithe & Turner 2002), particularly by stellar-size objects 
(microlensing). While the systematic errors introduced by K-correction, selection effects, and 
possible evolution can be removed, lensing might ultimately limit the accuracy of luminosity 
distance measurements (Perlmutter & Schmidt 2003). Only a large sample of SNe at each 
redshift can be used to characterize the lensing distribution and to correct for the effect of 
lensing. 

Properties of microlensed supernovae have been studied previously. Schneider & Wag- 
oner (1987) presented the time-dependent amplification of supernovae caused by the expan- 
sion of the photosphere, and showed that the related polarization of a supernova is not likely 
to exceed 1%. For cosmologically distant supernovae a 1% effect is practically impossible to 
detect. Linder, Schneider, & Wagoner (1988) studied amplification of supernovae and de- 
veloped approximate formulae for the amplification probability distribution. Rauch (1991) 
studied microlensing of SNe la by compact objects and calculated the resulting amplification 
probability distributions using Monte Carlo simulations. Kolatt & Bartelmann (1998) calcu- 
lated the light curves of type la supernovae microlensed by intracluster MACHO's assuming 
the point-deflector (with shear) model of Chang & Refsdal (1984). 

In this paper, we demonstrate how microlensing by a single stellar-size deflector can 
affect light curves of cosmologically distant SNe la. We use the simple Schwarzschild lens 
model to calculate amplifications. We ignore the additional amplification caused by the 
macrolensing introduced by the hosting galaxy, and concentrate on the time dependent 
effects caused by a single moving stellar-size deflector. We model the SNe la by expanding 
light sources with realistic radial surface intensity profiles. 

In the next section, we present a model light curve for type la supernovae. § 3 is 
devoted to a brief review of the gravitational lensing effect, and in § 4 we present the results 
of our calculations. Throughout, we assume a flat Friedmann-Lamaitre- Robertson- Walker 
cosmological model with Q m = 0.3, Qa = 0.7, and hi o = 0.67 to calculate the distances to 
the source and deflector. 
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2. Light Curves of Type la Supernovae 

In general, spectacular supernovae explosions can be interpreted in terms of two concepts 
(Arnett 1996): A shock wave running through a stellar envelope, and the radioactive decay 
of newly synthesized 56 Ni to 56 Co and then to 56 Fe. The nature of the explosion is the 
expansion away from the region of energy release. The explosion becomes more spherically 
symmetric with time, trending toward a small-scale Hubble type flow. 

In the meantime, various aspects of the event including shock emergence, radiative 
diffusion, heating, and recombination determine the evolution of the luminosity, i.e., the 
light curve, which is usually exhibited as the absolute magnitude as a function of time. In 
the case of SNe la, it is the shape of the light curve that can be used to measure cosmological 
parameters such as Hubble constant; see, for instance, Riess, Press, & Kirshner (1995). 

To model the intrinsic Type la supernova light curve we used a combination of two 
existing analytical models; for the peak of the light curve (early time) we used the model 
of Arnett (1982) while for the tail of the light curve (late time) we used the deposition 
function described by Jeffery (1999). The basic assumptions made in this model are: a 
homologous expansion of the photosphere, radiation pressure dominates at early times, the 
diffusion approximation is valid at early times, the optical opacity is constant for the light 
curve peak and the gamma-ray opacity is constant for the tail, 56 Ni is present in the ejecta 
and its distribution is peaked toward the center of the ejected mass at small initial radius. 
The details of this combined model are in Richardson, Branch, & Baron (2005). The model 
parameters have been fixed for a typical SN la: the kinetic energy is 1 foe (10 51 erg), the total 
ejected mass is 1.4 M , and the 56 Ni mass has been set to 0.6 M . These values produce a 
SN la with a peak absolute magnitude of about -19.5. 

To test the model we varied the kinetic energy, nickel mass, and the rise time until we 
achieved the best x 2 fit to the observed data of a few SNe la. The model worked well as 
can be seen in Figure 1, where the fit of SN1990N is shown with E k = 0.60 foe, M m = 0.62 
M Q , t rise = 21 days, and the reduced \ 2 — 0-43. We use this supernova because it was well 
observed (the light curve had good coverage) and it is photometrically characteristic of most 
SNe la. 



-4- 



3. Microlensing of Extended Sources 
3.1. Basics 

The linearized Einstein theory for a static gravitational field gives a bending angle for 
light rays passing through a weak gravitational field of 

a = — / V(f)dt, (1) 

J — oo 

where is the Newtonian gravitational potential satisfying the boundary conditions 0^0 
at infinity, and where the integral is performed along the light path in the absence of the 
gravitational field. For a point deflector the bending angle is simply 

AGm d . / v 

ol = — r, 2 

err 

where r is the light ray's impact vector and is the deflector's mass. Bourassa, Kantowski, 
& Norton (1973), and Bourassa & Kantowski (1974, 1976) used the 2-component nature of 
a to replace it with the complex scattering function I(z), where z = x + iy is the complex 
equivalent of the 2-d vector r = xi + yj. Using I(z), the relation between the source position, 
z, and image position, z , when both are projected onto the plane of the deflector (also called 
the sky plane) is 

z = z Q —I (z ) . (3) 

<r 

The scaled (effective) distance D is defined as D = Dd s Dd/D s , where the deflector-source 
distance, D^s, the observer-deflector distance, D^, and the observer-source distance, D s , are 
all the same type distances, e.g., apparent size distances. 

In the case of a point deflector (on which we focus in this paper) , the scattering function 
takes a very simple form: 

I{*o) = -, (4) 

Zo 

and the equation (3) reduces to 

z = z - -j, (5) 

z* 

where r E = \fAGm d D c~ 2 = ^/2r s D is the Einstein ring radius, and r s is the deflector's 
Schwarzschild radius. This equation has two separate solutions (images) for any source 
position r = \z\, 



r P= \zp\ = 2 ( \/^ 2 + 4r| + r 
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and 

1 



, , ., :1 ■ 14 '• ) • (6) 

Both images are in line with source and deflector. The 'primary' image, r p , lies on the same 
side of the deflector while the 'secondary' image, r s , is on the other side (Fig. 2). In the case 
of microlensing, the angular separation of the two images (9 P + 9 S in Fig. 2) is of the order 
of micro arcseconds. Consequently, the two images are seen as a single object (§ 3.2). The 
Einstein ring occurs when source and deflector are aligned with the observer (r = \z\ = 0) 
for which r p = r s = r E , and due to the symmetry of the lensing configuration the image 
is actually a ring. For a point deflector, the primary image is outside the Einstein ring 
ij'p ^ t e) while the secondary image falls inside (r s ^ r E ) as can be seen in Figure 3. 

The effect of gravitational lensing on the apparent brightness of a distant source can be 
computed in various ways. For extended sources like supernovae it is best to employ the fact 
that the apparent brightness is proportional to the image's apparent area, meaning that the 
brightness of a source is amplified by a factor A: 

A=^, (7) 

where A is the area of the image and A is the area of the source both projected on the 
deflector plane. For a point mass deflector and a small point-like source the amplification of 
the primary and secondary images respectively can be written as 

r 2 + 2r| 1 



and 



2rv / r 2 + 4r | 2 



= r 2 + 2r| _ 1 
S 2rVr 2 + 4r| 2 " 1 ' 



If the images are unresolvable the combined amplification becomes 



A P + A S = r y 2 t 2r f 9 - 0) 



'\/r 2 + 4r 2 



E 



If the source is not small, i.e., if r/r E varies significantly across the source, differential 
amplification must be taken into account. Measuring this amplification from a single obser- 
vation is not possible since it is practically impossible to figure out the original source flux. 
However, if the luminosity of the source varies with time in a predictable way as with SNe or 
if the source is of constant brightness and the lens is moving with respect to the line of sight 
to the source [as with observation of bulge stars; see, for instance, Sumi et al. (2004)], the 
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amplification will change with time in a predictable way and it is possible to determine the 
amplification. Figure 4 shows amplification curves of a point source lensed by a spherically 
symmetric deflector with mass m d = 1 M Q . 

The relative time delay of the two images (which is due to both the geometrical path 
difference of the two light signals and the difference in the gravitational potential through 
which the light rays travel (Cooke & Kantowski 1975) turns out to be only a small fraction 
of a second as expected for microlensing. Such small time delays between the two images 
can therefore have no effect on microlensed supernova light curve when features change on 
time scales of days or greater. 



3.2. Amplification of an Extended Source 

As implied in the last section, to obtain the total flux received from an extended source, 
an integral of intensity across the source may be required: 

A _ J image ° / 1f .x 

A ~J TdA- m 

J source 

If the surface brightness is constant across the source or if there is no differential amplifi- 
cation, the net amplification is simply given by equation (7), where A can be the primary, 
secondary, or total image area, giving respectively the primary, secondary, or total amplifi- 
cation. In the case of a circular extended source with a projected radius of a at a distance 
/ from the center of a spherically symmetric deflector (a and / are measured in the deflector 
plane), the total area of the combined and unresolvable images is 



A<***> = [ 2 a(a + lam<p)Jl + - . d<p (11) 

} _tl y I + a 1 + 2al sin ip 

After a rather long calculation, the total amplification is found to be 

A[a, /, r E \ = 77 {^K (k) + fi 2 E (k) + fi 3 U (n, k)} , (12) 

(Witt & Mao 1994; Mao & Witt 1998) where K, E, and II are respectively the first, second, 
and third complete elliptic integral with 

16a/r| 

k = 



(l + a) 2 ((l-af + 4r%) ' 
4a/ 

n = ^ , 

(a + 2 
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and constants 



7/ 



27ra 2 y / (/-a) 2 + 4r| 



^ = (I - a) (a 2 -I 2 - 8r|) , 
A* 2 = (I + a) ((I - a) 2 + 4r|) , 

_ 4(l-a) (a 2 + r|) 

/^3 r— — • 

t + a 

Using equations (8) and (9), equation (12)gives amplifications for both the primary and 
secondary images. Figure 5 is a plot of the amplification curves of an extended source 
moving with speed 0.5 AU per day (1 AU day -1 ~ 1.73 x 10 3 km s _1 ) with respect to a 
deflector of mass 1 M . To illustrate combined effect of relative motion and finite size on 
amplification, in Figure 5 we have assumed a constant surface brightness across the source 
(no limb darkening) and a fixed radius of 178 AU, corresponding to that of a supernova with 
a maximum atmospheric speed of 30,000 kilometer per second at eighteen days after the 
explosion. In Figure 5, the source approaches the deflector with various impact parameters 
(closest approach), producing different amplification curves. 

In what follows we use a more realistic supernova model to calculate amplification curves. 
We assume a ring-like structures for the intensity / of a spherically expanding photosphere 
like that of a SN la (Hoflich 1990). To model the 2-dimensional brightness profile across 
the source we assumed that the emissivity of light coming from any volume element in the 
photosphere of a type la supernova is constant. We found 



. 2 

r 



'M = Wi-^-j . (14) 

where I Q is the intensity at the center and r ph is the radius of the photosphere. This profile is 
in agreement with one of the commonly used family of limb-darkening profiles (Allen 1973; 
Claret, Diaz-Coroves, & Gimemez 1995). The assumption that emissivity per volume in SNe 
la is constant is a reasonable assumption during the nebular phase (t > 150 days). To see 
whether the obtained intensity profile is also applicable to photospheric phase (t < 150 days) 
we compared it to a normalized intensity distribution function (IDF) calculated using W7 
model (Nomoto, Thielemann, & Yokoi 1984) for SNe type la in U and B bands (courtesy of 
E. Lentz) . The result of this IDF calculation as well as our intensity profile can be seen in 
Figure 6. The two models show a reasonable match. 

We also estimated the characteristic radius of the optical image of a type la supernova 
as a function of time. For the photospheric phase, we assumed homologous expansion (r=vt) 
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and obtained the radius of the photosphere as a function of time since explosion by multi- 
plying the speed at the photosphere, as determined empirically by Branch et al. (2005), by 
the time. For the nebular phase, we assumed that the radius of the iron-group core expands 
at a constant velocity of 6000 km s -1 . A good fit for the speed is an exponential: 

u = 91e (-t/36.5days ) + ao? (15) 

in units of 10 3 km s -1 . Figure 7 shows the expansion velocity as well as the radius as a 
function of time. 



3.3. Probability 

The relevant quantity in seeing a microlensing event is the optical depth r. It is defined 
as the probability that a point source (or equivalently the center of an extended source) 
falls inside the Einstein ring of some deflector. The brightness of a point source within 
is amplified by a factor of at least 1.34. For randomly located point deflectors the optical 
depth depends on the mass density of the deflectors and not on their number density (Press 
& Gunn 1973). 



Typical values of the optical depths for microlensing of nearby luminous sources are re- 
markably small. For instance, Sumi et al. (2003) calculated an optical depth r = 2. 59^54 x 
10~ 6 toward the Galactic Bulge (GB) in Baade's window for events with time scales between 
0.3 and 200 days. Because of the small value of r, millions of stars need to be monitored 
when searching for microlensing in areas such as the GB, the Large Magellanic Cloud, or the 
Small Magellanic Cloud. The value of r is much higher for cosmologically distant sources. 
Assuming that the ordinary stellar populations of galaxies are the dominant causes of mi- 
crolensing events, Wyithe & Turner (2002) concluded that in a flat universe, at least 1% 
of high-redshift sources (z s ^ 1) are microlensed by stars at any given time. Zakharov, 
Popovic, & Jovanovic (2004) estimated that the optical depth for microlensing caused by 
deflectors both localized in galaxies and distributed uniformly, might reach 10% for sources 
at z s ~ 2. If we look at bulge-bulge lensing , Han & Gould (2003) give a model for the bulge 
from which a value of r = 0.98 x 10~ 6 is computed. They additionally compute an effective 
column density for deflectors = 2086 M pc~ 2 and a characteristic source-lens separation 
D = 782 pc defined by 

r = ^KD. (16) 

cr 

If we simply look at the bulge of a similar galaxy at z = 0.05 we expect a similar but 
D becomes the distance to the deflector D d and hence increased by a factor 2.7 x 10 5 . 
This would bring the optical depth up to 27% when looking through such a galaxy. At 
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this distance the size of the bulge is ~ 1 arcsec and clearly resolvable. The downside is 
one of alignment. What is the chance of a host galaxy being appropriately aligned with a 
foreground galaxy? The best place to see this effect seems to be the foregrounds of dense 
clusters. Because the typical galaxy is expected to contain a SNIa every 10 3 yrs, some 
3700(1 + z) alignments would have to be followed for a year to see one event. 

As mentioned above, the amplification of a point source falling inside the Einstein ring 
te is larger than 1.34. The probability of a larger amplification is proportionally smaller. 
The probability of having an amplification larger than A for a given lensing configuration is 

p(A) = u 2 A r(z s ) , (17) 

where ua = &A/ r s is the normalized impact parameter, resulting in amplification A (Paczyhski 
1986a,b) and t(z s ) is the optical depth for a source at redshift z s . For a point source, the 
result is 

u ^Htw=i- 1 )' (18) 

In the next section we use (17) to reduce the optical depth to probabilities appropriate for 
higher amplifications. 



4. Microlensed Light Curves 

In this section we use the model light curve of § 2 and the microlensing theory of § 
3 to predict the shape of lensed light curves of type la supernovae. We have calculated 
absolute magnitudes of the lensed SNe type la in the V-band. We assume that the source is 
at redshift z s = 1.0, hence introducing a time dilation factor of (z + 1) _1 = 0.5, and is lensed 
by a deflector located at redshift Zd = 0.05. Figure 8 shows the geometry of our model. To 
calculate the amplification as a function of time (eq. [12]), we need to know the distance / 
between the supernova and the deflector, projected on the plane of the deflector, 

lit) = yjll + vt («t T 2v^) , (19) 

where the minus sign is used when the supernova source explodes (t — and I = l Q ) before 
getting to the point of closest approach, 6, and the plus sign when it explodes after. We also 
need the Einstein ring radius te which is determined by the mass of the deflector and 
the distances of the supernova and the deflector from the observer. 

We plot light curves for various values of parameters b, l a , and v (the relative speed 
of source and deflector projected on the deflector's plane). Our sample deflector masses are 
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10 3 M Q , 1 M Q , and 10 M . We take the source to move in the deflector plane for 1,000 
days (observer time) with two relative projected speeds of the source v, 0.1 AU day -1 and 
0.5 AU day" 1 . 

Figures 9 through 17 show the amplification curves (upper panels) and the light curves 
(lower panels) of lensed supernova for nine different sets of input parameters. These inter- 
esting cases have been selected from a large number of configurations. The top panel of each 
figure shows the amplification of the primary and secondary images as well as their sum 
(total amplification), and the bottom panel includes light curves of the two images and their 
sum (the apparent light curve) as well as the original (unlensed) light curve of the supernova. 

We are plotting amplified absolute magnitudes of the supernova in the V-band, however, 
due to the redshift of the source, these light curves would be observed in the I-band. With 
the source located at z s — 1, My > —17 is too dim to be seen. Nonetheless, we include 
the complete amplification and light curves to show their trends over a period of 1,000 days 
after the supernova explosion. 

We can easily match the features in the light curve in each figure to the corresponding 
peak of the amplification curve. Figure 9 shows a huge overall increase in the brightness of 
the whole light curve. For this case the SN explosion occurs when the SN is near minimum 
impact l — b — 0. The amplification is strongest around the peak so the largest decrease 
in My happens there (almost 5 magnitudes). If the mass of the deflector were increased to 
10 M the magnitude of the 'narrowed' peak would reach My ~ —26 which is comparable 
to the magnitude of a bright quasar. Such huge amplifications can cause a bias in observing 
supernovae, allowing one to observe more distant supernovae, and as a result, to increase 
the depth (volume) of any supernova survey. Notice that this light curve looks different 
from those of strongly lensed supernovae in which, the lensed curve merely has an overall 
magnitude shift upward due to the static magnification introduced by the microlens' parent 
galaxy. 

In Figures 10 and 11, the lensed light curve hits a plateau before falling back on the 
expected trend of a type la supernova. For some values of the input parameters it is possible 
to observe a second peak in the light curve, e.g., see Figures 12 through 17. The second peak 
in each of these figures is bright enough to be observed and could easily be distinguished 
from the original (first) peak. If the impact parameter in Figure 12 is reduced to b = the 
second peak would actually be larger than the first. Figure 17 is another example of a second 
peak in the light curve which is bright enough to be seen although it occurs more than 400 
days after the supernova explosion. If the position of the initial explosion were increased to 
lo = 0.1 the still visible second peak would occur 2 years after the first. For many cases, 
microlensing has a less dramatic effect on the light curve's shape; it simply provides an 
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overall increase in its magnitude, and would be difficult to distinguish from amplification 
due to the galaxy hosting the microlens. 

The interesting cases described above occurred because the supernova exploded within 
the Einstein ring of a deflector, while the deflector proceeded to move on a time scale com- 
parable to the life of the supernova. Schneider & Wagoner (1987) found, in some instances, 
similar double peak modifications to the light curve caused by the supernova's photosphere 
expanding into a deflector's critical point. 

The standard optical depth r significantly overestimates the probability that the above 
interesting cases will occur. Using § 3 we can correct for the overestimate. For these lensing 
configurations, the normalized impact parameter ua does not exceed 0.1 (A ps 9) which, 
using equation (17), gives a maximum probability of ~ 10~ 4 for z s is 1 (r = 0.01), and 
~ 10 -3 for z s = 2 (r = 0.1) if Zakharov, Popovic, & Jovanovic (2004) are correct. These 
numbers are somewhat higher than Linder, Schneider, & Wagoner (1988) were predicting 
for Type 1 SN but not for Type II. With a supernova rate of R SN eia = 2, 000 yr" 1 the SN 
lensing rate at redshift z = 1 is ~ 10 4 S1 yr _1 (Oguri, Suto, & Turner 2003). Rates likes 
this, together with time scales of some cases studied here, imply that such effects may not 
be observed unless a large number of cosmo logically distant (around 10 4 for z s ^ 1) type la 
supernovae are followed for a period of up to 2 years. For the high amplification cases with 
a second peak, a seperate probability estimate can be made (see the Appendix). 



5. Conclusion 

We have shown that microlensing can significantly affect the light curve of a cosmologi- 
cally distant type la supernova. We restricted our calculation to z s — 1 and Zd = 0.05 in the 
currently accepted Vt m = 0.3, Q\ = 0.7 flat cosmological model. We found that microlens- 
ing not only increases the magnitude of the light curve but also can cause a change in its 
shape. Relative transverse motion of the SN and lens, as well as the expanding photosphere 
(Schneider & Wagoner 1987) can result in features such as a narrow but high peak, a plateau 
following the peak, or even the presence of a second peak. 

In § 4 we concluded that for microlensing by compact masses distributed through the 
cosmos, the optical depth is 0.01 (z s ~ 1) but might reach 0.1 (z s ~ 2), meaning that the 
overall chance of a distant supernova type la being microlensed is not negligible. Any multi- 
band supernova survey aimed at finding supernovae at redshifts around z — 1 (and above), 
could discover and identify one microlensed SN la event out of roughly a hundred events. 
However, the low impact parameters required to produce the special features depicted in § 
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4 demand observation of ~ 10 4 supernovae at z s ^ 1. And, to see unusual features such as 
double peaks, the lensed supernova must be followed for an extended period of ~ 2 years. In 
the Appendix we have made an optical depth type estimate to include double peak events. 
For microlensing by stars in the bulge of a galaxy at z d = 0.05 we find a max probability of 
~ 1.7 x 10~ 3 . This is an ideal deflector distance for observing double peaks due to transverse 
motion. As expected this number is only slightly smaller than the 27% optical depth estimate 
made in §3.3 for bulge lensing when corrected for an impact u = 0.1 These estimates do not 
take into account the observational bias in favor of amplified events (§ 4) nor the possibly 
enhanced probability due to evolution of the rate at which Type-la supernovae have occurred. 
It is interesting to note that according to OGLE III (Udalski 2003) a histogram of u values 
for Bulge lensing peaks at u ~ 0.1 (probably due to amplification biasing). 

The authors are pleased to thank an anonymous referee for helpful criticisms of the 
first version of the manuscript and Zach Blankenship for correcting an oversight in this 
version. This work was in part supported by NSF grant AST0204771 and NASA grant 
NNG04GD36G. 



A. Appendix 

In this appendix we compute the probability that a source, followed for a period T, 
impacts a point deflector with a reduced impact parameter less than u = b/R,E and simul- 
taneously moves at least a distance b during the period T. Such a time dependent impact 
will cause a change in the amplification of 10%-50% depending on the actual impact. The 
idea here is to estimate the chance of seeing a distortion in the light curve of a SN whose 
life time is T$n ~ 200 days. 

We start with a number density N d of mass m deflectors (located at a distance D d from 
the observer) moving with relative transverse velocities distributed according to: 

f/Vd =A,™^e-^ ms . (Al) 



av v, 



2 

rms 



The probability of one of these moving deflectors impacting the line of sight to a source at D s 
with a reduced impact parameter ^ u and moving a reduced distance ^ u during a period 
T is: 



V _„,2 /o,,2 



AProb(u,T)= dD d I dvN d (D d )(nu 2 R 2 E + 2uR E vT)^—e- v2/2v r™. (A2) 

JO JuR E /T V rn 



The integrand is the sum over areas represented in Figure 18. The velocity integral can be 
done easily, and if the deflectors are effectively confined to a plane, the result can be written 

as 



AProb(u, = ^ £ D u 2 | (tt + 2)e^ 2 / 2 + J , (A3) 



where 

z = —T = u -r> (A4) 

u rms - 1 ± 

E is the projected surface mass density 

mN d (D d )dD d , (A5) 

./o 

and .EV/c is an error function. 



Jo 



The characteristic crossing time for microlensing is defined by T rms = RE/v rms = 
\/2rsD/v rms (see § 3 for definitions) which for Galaxy bulge-bulge lensing is about 10 days 
(Udalski 2003). For a similar galaxy at redshift z = 0.05 lensing a distant SN through its 
bulge, the reduced distance D is increased by a factor of ~ 2.7 x 10 5 [see § 3.3 and Han & 
Gould (2003)] and hence T rms increases to ~ 5,200 days. If u ~ 0.05 and T = T$n ~ 200 
days, then £ ~ 1.3 and APro6(0.05, 1.3) ~ 1.7 x 10 -3 . This particular probability falls off 
by at least an order of magnitude when u < 0.01 or u > 0.15. 
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Fig. 1. 

Model light curve of a supernova type la in its rest frame. Observed data of SN 1990N 
are added to show the fit. 

Fig. 2.- 

Microlensing by a point deflector. The deflector's position is considered as the coordinate 
origin in the deflector plane. 

Fig. 3.- 

The position of the primary and secondary images on the deflector plane (angular coor- 
dinates) with respect to source and Einstein ring. Note that source makes an angle of 0.04 
micro arcsecond with respect to the deflector. 

Fig. 4.- 

Amplification curves of images of an extended source (z s = 1.0) by a deflector at redshift 
z d = 0.5. The horizontal axis is the angular separation of deflector and source, (3, (see Fig. 
3) in micro arcsecond. 

Fig. 5. — 

Amplification curves of an extended source with a radius of 178 AU, moving at 0.5 AU 
day -1 in the deflector plane. The source reaches its closest approach (impact parameter) at 
t = 0. The figure shows five trajectories with different impact parameters, b (see Fig. 8). 
For this figure only, the surface brightness and radius are assumed constant. 

Fig. 6.— 

Normalized IDF profile calculated for SN la in U and B bands using W7 model (courtesy 
of Eric Lentz, University of Georgia), used to calculate the amplification curves of a supernova 
as an extended source with radius- dependent surface brightness. In this figure, a maximum 
expansion velocity of 30,000 km s _1 has been used to obtain the expansion velocity projected 
on the sky. The intensity profile used in the paper is also presented for photospheric speeds 
of 13,000 km s- 1 and 15,000 km s _1 . 

Fig. 7.- 

Expansion speed (upper panel) and the radius (lower panel) of the photosphere of a 
type la supernova as a function of time. 



Fig. 8.- 
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Geometry of the microlensing model used in § 4. 
Fig. 9- 

Alification curve (upper panel) and light curve (lower panel) of a supernova at z s = 1.0, 
microlensed by a deflector at z d = 0.05. The input parameters are m d = 1 M , v = 0.5 AU 
day -1 , 6 = 0.0 r E , and Z = 0.0 r^. 

Fig. 10.- 

Same as Fig. 9, with input parameters m d = 10 -3 M , v = 0.1 AU day -1 , b = 0.1 r E , 
and Z G = 0.25 r E . 

Fig. 11.— 

Same as Fig. 9, with input parameters m d = 1 M , v = 0.5 AU day -1 , b = 0.01 r E , and 
l = 0.05 r B . 

Fig. 12.- 

Same as Fig. 9, with input parameters = 10 -3 M Q , v = 0.1 AU day -1 , b = 0.1 rg, 
and l = 1.0 r E - 

Fig. 13.- 

Same as Fig. 9, with input parameters nid = 10 -3 M , v = 0.5 AU day -1 , b = 0.01 r E , 
and l = 1.5 r E . 

Fig. 14.— 

Same as Fig. 9, with input parameters = 10" 3 M , v = 0.5 AU day -1 , b = 0.1 r E , 
and l = 4.0 r E . 

Fig. 15.- 

Same as Fig. 9, with input parameters = 1 M , v = 0.5 AU day -1 , b = 0.0 r E , and 
l o = 0.1r E . 

Fig. 16.- 

Same as Fig. 9, with input parameters m d — l M , f = 0.5 AU day -1 , b = 0.01 r E , and 
Z o = 0.1r E . 

Fig. 17.- 

Same as Fig. 9, with input parameters m d = 10 M , v = 0.5 AU day -1 , b = 0.0 r E) and 
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l = 0.05 r E . 
Fig. 18.— 

Area on the deflector plane within which a microlens would be close enough to a luminous 
SNe la to cause significant changes in its lightcurve. 
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1 1 1 1 1 

Primarv Tmatrp 

Finstpin Ring 
Snnrrp 


| i | i | i | 

deflector coordinate: (0.0, 0.0) 

source coordinate: (0.0, 0.04) 

,.•* 


Secondary Image 
1,1, 


o _ 

1,1,1,1 



-3 



-2 -1 

angular position (\i arcsec) 



Figure 3 



100 


1 1 

i\ 

- '1 
i\ 
l\ 
i\ 
i\ 


1 1 1 1 1 1 1 

- - ■ primary amplification _ 
■ ■ ■ ■ secondary amplification z _ 

— total amplification : 


10 

1 

0.1 


i \ 
\ \ 

— \ \ 
•\ \ 
'.\ \ 




1 , 1 


I , I , I , I 



J I I I I I I I I I l_ 

0.5 1 1.5 2 2.5 



P ((0. arcsec) 



Figure 4 




Figure 5 




Figure 6 



-25 - 




t (days) 



Figure 7 



- 26 - 





1 Deflector ^ 


r 

E 




\ b \ / 






\ ri )\ , 




r : Deflector's Einstein 

E 


Ring Source 




1 : Distance Between Moving Source & Deflector 




1 : Distance Between Source & Deflector when SN explodes 




b : Impact Parameter 






V : Relative Transverse Speed of Source 




Figure 8 



-27- 




Figure 9 



-28- 




Figure 10 



-29- 




Figure 11 




Figure 12 



-31 - 




200 400 600 800 1000 

t (days) 



Figure 13 




Figure 14 



T 




200 400 600 800 1000 

t (days) 



Figure 15 



-34- 




Figure 16 



-35 - 



400 



300 - 



o 



I 200 



100 


-24 



-21 ' 
18 [f 



■15 - 



•12 
-9 



Total 

Primary 

Secondary 




V 



Apparent 
Primary 
Secondary 
Unlensed 



m = 10 M 

d sun 



- v = 0.5AUday 
. b = 0.0 r E 
1 =0.05r c 

o E 



1 







200 400 600 

t (days) 



800 



1000 



Figure 17 




Figure 18 



